[V- 

CL, 


o 

C_2 


Lefschetz  Center  for  Dynamical  Systems  ^ 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  Rhode  Island  02912 


A " 

I " ' 


(T/^Dec 


78 


/ 


/ j K!  1V-74  d. 


- — — \ 

1 ^ - Tol^  

\ tor  P»bU' 1 Vs  — 

1 ill  |i-;1-,lhoa  — 


This  research  was  supported  in  part  by  the  National  Science  Foundation 

in  Part  by  tbe  U.S.  Army  under  contract 
ARO- DAAG27 - 76-G-0294  and  An  part  by  the  Office  of  Naval  Research  under 
contract  N0014 - 76-C- 0278 -P0002 . 

G 


Vo±  0 1 


J 


1 


1.  Introduction 


^ There  are  several  indications  pointing  to  an  intimate  re- 
lationship between  the  second  law  of  thermodynamics  and  "stability". 
Most  notably,  the  work  of  Ericksen  [1]  on  thermoelasticity  and 
subsequent  developments  by  Coleman  and  Dill  [2]  and  Gurtin  [3] 
have  demonstrated  that  the  Clausius-Duhem  inequality  induces 
Liapunov  stability  of  equilibrium  processes.  ^Herc  I attempt  to 
establish  a different  connection  between  stability  and  the  second  law, 
in  that  "stability"  will  be  interpreted  as  continuous  dependence 
of  thermodynamic  processes  upon  initial  state  and  supply  terms. 

The  ideas  will  be  presented  within  the  context  of  thermo- 
elasticity theory,  without  heat  conduction.  The  hyperbolic  character 
of  the  field  equations  causes  the  breakdown  of  smooth  solutions  and 
the  development  of  shock  waves  so  the  class  of  smooth  functions  is  far 
too  narrow  to  encompass  all  processes  of  physical  interest.  The 
natural  framework  is  the  class  of  proper  processes,  characterized 
by  the  property  that  velocity,  deformation  gradient,  specific  entropy, 
stfess,  temperature  and  internal  energy  are  functions  of  bounded 
variation. an  the  sense  of  Tonelli-Cesari . J The  balance  laws, 
in  integral  as  well  as  in  local  form,  are  meaningful  in  the  class 
of  proper  processes  and  the  classical  geometric  theory  of  wave 
propagation  can  be  transplanted  into  this  class.  One  expects  that 


that  is,  measurable  functions  whose  first  derivatives  are  Borel 
measures.  This  is  essentially  the  broadest  function  class  in 
which  the  Gauss-Green  theorem  holds.  There  is  an  analogy  between 
the  geometric  structures  of  functions  of  bounded  variation  and 
functions  that  are  piecewise  smooth.  For  a survey  see  [4] . 
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existence  of  solutions  to  the  field  equations  will  eventually  be 
established  within  the  class  of  proper  processes  but  so  far  this 
has  only  been  accomplished  in  the  one-dimensional  case  [ 5 1 . 

The  interpretation  of  the  role  of  the  second  law  of  thermo- 
dynamics in  the  theory  of  proper  processes  poses  an  interesting 
problem.  The  idea  of  restricting  constitutive  relations  so  that 
the  Clausius-Duhem  inequality  be  automatically  satisfied  by  smooth 
processes  originated  in  the  work  of  Coleman  and  Noll  t6 1 and  has 
by  now  become  standard  practice  in  continuum  thermomechanics.  On 
the  other  hand,  within  the  class  of  processes  with  shock  waves,  the 
Clausius-Duhem  inequality  has  traditionally  been  viewed  as  an  ad- 
missibility criterion.  It  turns  out  that  this  criterion  generally 
rules  out  some  but  not  all  extraneous  processes  so  that  one  must  either 
impose  chock  admissibility  restrictions  [7]  or  else  strengthen  the 
second  law  (8 1 , in  order  to  single  out  those  processes  that  are 
physically  relevant. 

Despite  the  above  remarks,  it  will  be  shown  here  that,  when- 
ever they  exist,  smooth  processes  are  stable  within  the  class  of 
(not  necessarily  smooth)  proper  processes  that  satisfy  the  Clausius- 
Duhem  inequality.  In  other  words,  for  as  long  as  one  is  dealing 
with  smooth  processes,  the  second  law  in  its  traditional  form  is  a 
satisfactory  admissibility  criterion.  * 

The  above  result  is  established  under  certain  assumptions 
on  material  response,  relating  to  convexity  of  internal  energy. 

With  the  exception  of  one-dimensional  bodies,  global  convexity  of 
the  internal  energy  function  is  contrary  to  experience  (it  is 
incompatible,  in  particular,  with  the  principle  of  frame  indifference). 
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Even  so,  internal  energy  will  in  general  be  locally  convex  on  a 
certain  region  in  state  space.  In  Section  4 we  prove  uniqueness  and 
stability  for  smooth  processes  residing  in  the  convexity  region  of 
internal  energy.  For  the  case  where  the  body  is  smooth  and  bounded  and 
the  motion  of  its  boundary  is  prescribed,  stability  is  established, 
in  Section  5,  under  the  weaker  assumption  that  the  smooth  process 
resides  in  the  strong  ellipticity  region. 

The  general  strategy  of  the  proof  was  inspired  by  the 
important  paper  [9]  of  DiPerna  on  uniqueness  of  solutions  to  the 
initial  value  problem  for  quasilinear  hyperbolic  systems.  The  heart 
of  the  proof  is  an  inequality,  derived  in  Section  3,  which  estimates 
the  evolution  in  time  of  the  "distance"  between  the  states  of  two 
processes  originating  at  neighboring  states.  The  time  rate  of 
increase  of  this  distance  can  be  controlled,  with  the  help  of  the 
Clausius-Duhem  inequality,  provided  that  at  least  one  of  the  processes 
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2.  Adiabatic  Processes  in  Thermoelastic  Materials 


We  consider  a thermoelastic  body  with  reference 
configuration  £ in  Rn  (n  = 1,2  or  3).  The  reference  mass 
density  P(X)>  defined  on  is  smooth  and  strictly  positive, 

(2.1)  P (X)  > PQ  > 0,  X e ® . 

A motion  x = x(X,t)  determines  the  velocity  field  v = x and 
the  deformation  gradient  field  F = grad^x.  The  internal  energy  e, 
the  Piola-Kirchhof f stress  T and  the  temperature  0 are 
determined  by  the  deformation  gradient  F and  the  specific  entropy 
n via  constitutive  relations 


(2.2) 


e = e*(F,n;X),  t = T*(F,n;X) , 0 = 0*(f,p;X), 


where  e*,T* 
set  M+  of 
and  X in 


•j^ 

and  0 are  smooth  functions  defined  for  F in  the 
n x n matrices  with  positive  determinant,  p in  R 

75-  (*) 

& • We  also  require  that 


(2-3)  e*(F»n;X)  > 0,  F e M+,  n E R,  X c 

\ 

We  shall  be  assuming  that  the  material  is  a non-conductor  of  heat  so 
that  heat  flux  and  entropy  flux  vanish. 


pn ■ ■ — 

1 n particular,  we  assume  that  the  partial  derivatives  of  e*,0* 

and  T , at  any  fixed  F e M and  p c R,  are  bounded  functions 
of  X on  9. 
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By  a proper  thermodynamic  process  we  mean  fields  (x,n)(X,t) 
such  that  (v,F,n)(X,t)  are  functions  of  locally  bounded  variation, 
in  the  sense  of  Tonel 1 i -Cesar i , satisfying  the  balance  laws  of 
momentum  and  energy,  viz., 


(2.4) 


Pv. 

l 


T . + 

ia,a 


Pbi, 


(2.5)  i v.vp  - (Tiov.)iCl  * ebjVj  ♦ pr, 

where  b(X,t)  is  the  body  force  and  r(X,t)  is  the  energy  supply. 
The  process  will  be  called  admissible  if  it  also  satisfies  the 
Clausius-Duhem  inequality 


(2.6)  n - - > 0. 

0 ' 

The  reader  should  keep  in  mind  that  for  the  class  of  proper  processes 
(2.4)-(2.6)  only  hold  in  the  sense  of  measures  (or  distributions) 
and  thus  one  cannot  simplify  (2.5)  in  the  standard  fashion,  since 
the  usual  product  differentiation  rule  does  not  apply  to  products 
of  functions  of  bounded  variation. 

A process  (x,7})(X,t)  will  be  called  smooth  if  the  functions 
(v,F,  rf)  (X,t)  are  Lipschitz  continuous,  uniformly  on  bounded  subsets 
of  their  domain.  Thus  a smooth  process  may  contain  weak  waves  but 
not  shock  waves.  For  smooth  processes,  one  may  write  the  balance 
laws  in  reduced  form, 


(2.7) 


Pv  . = T . + pR . . 

vi  1 KDi» 


M - 


rmsm. 


As  it  is  known,  the  standard  requirement  that  every  smooth  process 
be  admissible  will  he  satisfied  if  and  only  if 


(2.9) 


dt7 

3F 


3e 

3n 


Indeed,  upon  using  (2.9),  (2.8)  yields 


(2.10) 


0, 


which  shows  that  the  Clausius-Duhem  inequality  is  automatically 
satisfied  by  smooth  processes,  as  an  equality. 
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3.  The  Evolutionary  Inequality 

We  establish  here  the  estimate  that  will  be  the  basis  of  our 
stability  analysis.  We  define 


(3.1)  H* (v,F, n;v,F, n;X)  = j P(vi-v.) (v.-v.)  + pe*(F,n;X) 


- pe*(F,n;X) 


po* (f, n;X) (n- n) , 


(3.2)  G;(v,F,ri;v,F,n;X)  = -(T?a(F,n;X)  - T*u(F,n;X))(v.-vi) 


On  account  of  (2.9),  H and  G are  of  quadratic  order  in 
(v-y,F-F,n-n) . 

Consider  now  a smooth  process  (x,n)(X,t),  with  supply  terms 
(b,r)(X,t),  and  let  (x,n)(X,t)  be  any  admissible  process  with 
supply  terms  (b,r)(X,t).  We  set 


(3.3)  H (X , t)  = H (v(X,t),F(X,t),n(X,t);v(X, t),F(X,t),n(X,t);X), 


(3.4)  G(X,t)  = G*(v(X,t) ,F(X,t) ,n(X,t) ;v(X,t) ,F(X,t) ,n(X,t) ; X ) . 


We  shall  be  viewing  H as  a measure  of  the  "distance^’  between  the 
two  processes.  In  order  to  see  how  H evolves  in  time,  we  compute 
below  H + div^G,  keeping  in  mind  that  the  usual  product  differentiation 
rule  applies  (in  the  sense  of  measures)  to  the  product  of  a Lipschitz 
continuous  function  with  a function  of  bounded  variation.  We  have 


.■mm. 


(3.5)  II  + Gu>ot  = P(e  + 2 ViVi)  ‘ PViVi 


PV . V . + pv . V . 
11  11 


* itp  " + *iu*ia  ’ 


- P6Cn-Tt)  • (T^vj)  * CTl0-Tl0)F10  * T.^ 


+ v . T . + v . T . - v . T . 

l ia,a  l i«,a  i lU.a1 


Observing  that  v.  a = F^a,  a = F\  a and  using  the  balance  laws 
(2.4),  (2.5),  (2.7),  (2.8),  we  may  rewrite  (3.5)  in  the  form 


(3.6)  H ♦ Ga>a  = P(v.-v.)(b.-Fi)  ♦ P(r-r)  - Tj3(Fjp-Fjp) 


+ Fia(Tia-Tia)  ' P0Cn-n)  - P0(n-n) 


P(vi-vi)(bi-bi)  - P(r-r)  - Fiu(F.g-Fjg) 


3T"  • _ i _ an*  .L 

n(Fje-Fj3)  + Fiadia-Tia)  - p |p-  Fia(n-n) 


30  j.  _ _ . jl 

p z7r  n(n-n)  - po (n-n) . - 


By  virtue  of  (2.9), 


(3.7) 


3Fia  - 3n  • ■ 3Fj6  • 


so  that  (3.6)  yields 


T~ 
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(3.8)  II  + Gu>ct  = P(v.-vi)(b.-hi)  + P(r-r) 


^ _ 3T*  _ 3T* 

* >VTi«  ' V - nr;  (Fjs'Fje)  ' jf 


. pTJte  - ir  - |^-  (Fje-rje)  - (n-n)) 


- pn(0-B’)  - p^(n-n) 


From  the  Clausius -Duhem  inequality  (2.6)  and  (2.10),  we  obtain 


(3.9)  P(r-r)  - Pn(9-0")  - pH'(n-n)  < £ (9-S')(r-r)  - — (« -IT) 2 

' a 00“ 


hence  (3.8)  takes  the  form 


(3.10) 


H + G < P(vi  -V:)  (b. -b.)  + £ (B-0")(r-r)  - — (9-0)2 
a,a  i 1 1 1 u 0tT 


3T*  _ 3T* 

ia  1 * ia  x ia  9FT^  (F  j 3_F  j 3^  ' JrT~  (n'n)} 


+ F.„(T.„  - T . . 


+ f<e  3-F3 


Inequality  (3.10)  will  be  the  starting  point  of  our  stability 
analysis.  The  crucial  observation  is  that  its  right  liand  side  is  of 
quadratic  order  in  (v- v , F- F , n- 0 , b-b , r- r)  so  that  one  may  establish 
continuous  dependence  of  processes  upon  initial  data  and  supply  terms 
by  applying  Gronwall  type  inequalities,  provided  that  H be  positive 
definite  in  an  appropriate  sense.  This  program  will  be  implemented 
in  the  following  two  sections. 


4.  Stability  of  Smooth  Processes  in  the  Convexity  Region 

of  Internal  Energy 
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In  this  section  we  shall  be  dealing  with  pairs  of  processes 
that  render  H(X,t)  pointwise  positive  definite.  On  account  of 
(3.1)  and  (2.9),  the  sign  of  H is  dependent  upon  the  convexity 

'ft 

properties  of  the  function  e . 

For  fixed  X e 3 , we  let  J*(F,n;X)  denote  the  Hessian 
matrix  of  e*(F,n;X).  We  will  say  that  a process  (x,n)(X,t) 
resides  in  the  convexity  region  of  internal  energy  if 

ft  

J (F(X, t) , n(X, t) ;X)  is  uniformly  positive  definite  on  the  closure 
of  the  domain  of  the  process.  By  virtue  of  (2.9),  J*  is  the 
Jacobian  matrix  of  the  transformation  (F,n)  H-  (P  ^T,9).  Thus, 

ft  __  — 

by  Sylvester's  theorem,  J (F (X , t) , n(X, t) ; X)  will  be  positive 
definite  if  and  only  if 


(4.1) 


3T*(F(X,t),n(X,t) ;X) 


9F 


positive  definite 


(4.2) 


r*  ,77 


det  J(F(X,t) ,n(X,t) ;X)  > 0. 


In  order  to  see  the  thermodynamic  interpretation  of  f4f.2),  we 
observe  that  when  J*  is  nonsingular  at  a point  one  may  invert 

locally  the  transformation  (F,n)  h-  (P  ^T,9)  to  get 

* - 1 * . 1 

i:  — ■Lrr  ~ _ r ^ x-r 


an 


. -n , 


(F, n)  K 

(P_1 

T, n;X) 

and 

f 

*)  ^det 

9I 

3F 

(4.3) 
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It  follows  that  (4.1),  (4.2)  are  equivalent  to  the  following  two 
conditions  (in  the  notation  of  classical  thermodynamics): 


(4.4) 


(4.5) 


9T 


9F 


aM 

3fJ. 


positive  definite, 


> 0, 


It  is  easily  seen  that  (4.4),  (4.5)  are  in  turn  equivalent  to 


(4.6) 


positive  definite, 


(4.7)  (l?)E  > »• 

The  relevance  of  (4.5),  (4.7)  to  thermodynamic  stability  was  first 
pointed  out  by  Gibbs  [10].  It  would  be  appropriate  to  impose  these 
conditions  on  any  region  in  state  space  that  does  not  contain 
critical  points  at  which  the  material  undergoes  phase  transitions. 

On  the  contrary,  (4.4)  and  (4.6)  are  quite  restrictive.  These 
conditions  may  be  globally  satisfied  for  one  dimensional  bodies  but 
are  incompatible, in  the  multidimensional  case,  with  the  principle 
of  material  frame  indifference  on  a large  portion  of  ^tate  space 
that  includes,  in  particular,  natural  states.  The  theory  of  elastic 
stability  (e.g.  [11])  elucidates  the  role  of  (4.4),  (4.6)  and  makes 
it  plausible  that  these  conditions  will  be  satisfied  on  some  region 
in  state  space  where  strains  are  moderate  and  stresses  are  pre- 
dominately tensile. 
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We  now  state  two  stability  results,  the  first  one  for  the 
mixed  initial-boundary  value  problem  and  the  second  for  the  pure 
initial  value  problem. 

Theorem  4.1.  Assume  that  0 is  bounded  and  has  finite  perimeter 
[4],  Let  (x,  n) (X , t ) be  a smooth  process  defined  on  & * [0,tp], 
residing  in  the  convexity  region  of  internal  energy,  with  supply 
terms  (b,F)(X,t)  e L°° ( & * [ 0 , t Q ] ) . Then  there  are  positive 
constants  6,oi,M,N  with  the  following  property:  If  (x,n)(X,t) 
is  any  admissible  process  defined  on  & x [ 0 , t ^ ] , with  supply  terms 
(b , r)  ( • , t)  e L1([0,tQ];  L2(.0)),  and  such  that*- 

(4.8)  | F (X, t)  - F(X,t) I + In(X.t)  - n(x,t)|  < 6,  (X,t)  e m x [0,t0l, 

(4.9)  (v-v).(T-T)  < 0,  on  8.0, 
then  we  have,  for  any  s e [0,tg], 

(4.10)  I I (v-v,F-F,n-n) ( - , s) | | 2 

L (.0) 

n c 

< Me  || (v-v,F-F,n-n) (• ,0) | | . 

L (Q  \ 

+ NeUS  fS| | (b-b,r-r) (• ,t) | | - dt. 

J0  LZ(0) 

T3Tj Z 

v -'In  Equation  (4.9)  t,t  denote  the  stress  vectors. 


I 

i 

i 
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Theorem  4.2.  Assume  that  & = Rn.  Let  (x,n)(X,t)  be  a smooth 
process  defined  on  & x [0,tQ],  residing  in  the  convexity  region 
of  internal  energy,  with  supply  terms  (b,r)(X,t)  e L^oc^  x [ 0 , t ^ ] ) 
and  such  that  (F,n)(X,t)  is  bounded  on  x [ 0 , t q ] . Then  there 

are  positive  constants  6,k,ot,M,N  with  the  following  property:  If 
(x,  n)  (X,  t)  is  any  admissible  process  defined  on  ^ x [ 0 , t ^ ] , with 
supply  terms  (b,r)(*,t)  e L^([0,tg];  L^oc(3#)),  and  such  that 


(4.11)  |F(X,t)  - F(X,t) | + |n(x,t)  - n(X,t) I < 6,  (X,t)  E ,^x  [0,tQ] , 


then  we  have  for  every  a > 0,  s e [0,tg], 


(4.12)  | | (v-v,F-F,n-n)(-,s)|  | 2 

LZ ( | X | <a) 


< Me  1 1 (v-v,F-F,n-n) (• ,0) 1 1 7 

( | X | <a+ks) 


QtS  fS 

Ne  | | (b-b,r-r) (• ,t) | | 2 dt. 

J0  ' L^(|X|<a+k(s-t)) 


From  the  above  propositions  one  draws  immediately  the  following 


Corollary  4.1.  Let  (x,"n)(X,t)  and  (x,n)(X,t)  be  as  in 

Tlicorem  4.1  or  as  in  Theorem  4.2.  Assume  that  the  corresponding 
supply  terms  (fr,r)(X,t)  and  (b,r)(X,t)  coincide  on  @ x [ 0 , t Q ] 
and  that  both  processes  originate  from  the  same  state,  that  is, 


(4.13) 


x (X , 0)  = x (X , 0)  , v(X,  0)  = v (X, 0)  ,.  X £ & 


I 
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Then  (x,n)(X,t)  and  (x,n)(X,t)  coincide  on  ^ * [ 0 , t ^ ] . 

The  above  uniqueness  result  is  only  local  since  the  two 
processes  are  restricted  a priori  by  (4.8)  or  (4.11).  Roughly 
speaking,  the  Corollary  states  that  the  only  way  that  an  admissible 
process  may  bifurcate  out  of  a smooth  process,  at  an  instant  t^, 
is  by  the  spontaneous  generation  at  t^  of  a shock  wave  of  large 
(>  6)  amplitude.  In  the  one  dimensional  case  it  is  known  [12]  that 
this  is  impossible,  i.e.,  shock  waves  originating  at  a smooth  state 
start  our  with  "zero"  amplitude  and  are  then  gradually  amplified. 

One  would  need  a similar  property  for  the  multidimensional  case  in 
order  to  deduce  global  uniqueness  from  Corollary  4.1  but  a proof 
does  not  seem  possible  at  the  present  stage  of  development  of  the 
theory. 

For  the  proof  of  Theorems  4.1  and  4 . 2 we  shall  need  the 
following  Gronwall  type  inequality. 

00 

Lemma  4.1.  Assume  that  the  nonnegative  functions  y(t)  £ L [0,s] 
and  g(t)  £ L 1 [ 0 , s ] satisfy  the  inequality 

(4.14)  y2(o)  < M2y2 (0)  + f°[2ay2(t)  + 2Mg , t) y (t) ) dt , ae  [0,s], 

0 

i 

where  u,m,N  are  nonnegative  constants.  Then 

(4.15)  y(s)  < MeU  y ( 0 ) + NeUS  g(t)dt. 

0 

The  proof  of  the  above  lemma  is  straightforward  and  in  any 
case  a more  general  result  (Lemma  5.1)  will  be  established  in  Section  5. 
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) 

1 

! 

► 

Proof  of  Theorem  4.1.  Since  (x,n)(X,t)  resides  in  the  convexity 
region  of  internal  energy,  there  exist,  on  account  of  (2.9), 
positive  constants  6, A such  that,  whenever  (4.8)  holds, 

(4.16)  e*(F(X,t),n(X,t);X)  - s*  (F  (X,  t)  , o(X,  t)  ;X) 

- P'1T*a(F(X,t),n(X,t)  ;X)[F.a(X,t)  - F.JX.t)] 

- 0*(F(x,t) ,n(x,t) ;x) [n(X,t)  - n(x,t)] 

> A(|F(x,t)  - F(X,t)|2  + | n(x,t)  - n(x, t) | 2) , 
for  every  (X,t)  e <2\  * [0,tQ). 

We  now  fix  s e [ 0 , t Q , we  integrate  (3.10)  over  & x 1 0 , a) , 
a e [0,s],  we  apply  the  Gauss-Green  theorem  and  we  use  (4.9),  (4.16) 
and  the  observation  that  the  right-hand  side  of  (3.10)  is  of 
quadratic  order  in  (v-v , F-F , n- ri, b-b , r-r)  thus  arriving  at  an 
estimate  of  the  form  (4.14)  with 


(4.17) 

y(t)  = 

| | (v-v,F-F,n-n) ( * , t)  | | ? 

L^(<3) 

(4.18) 

g(t)  = 

| | (b-F,r-r) (• ,t) | | 2 

LZ(#) 

An  application  of  Lemma  4.1  completes  the  proof. 


Proof  of  Theorem  4.2.  As  in  the  proof  of  Theorem  4.1  there  are 
positive  constants  6, A such  that  (4.16)  is  satisfied  for  all 
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(X,t)  e . <3f  * 10,  t]  , whenever  (4.11)  holds.  Consequently,  since 

G*  is  of  quadratic  order  in  (v- v , F-F,  n- fi)  , one  may  determine  a 
sufficiently  large  positive  constant  k with  the  property 

X 


(4. 

19) 

kH  ( X , 

t) 

- Q(X,t) 

■ m - 

0, 

(X,t)  e <2  * 

to, t0i 

• 

We 

now 

fi 

x s £ 

[0 

» f Q J * 3 

> 0,  and 

for 

each  a e [0 

, s ] we 

i ntegra tc 

(3. 

10) 

over  the 

f rust urn 

{(X,t)  | 

t c 

[0,o],  | X | < 

a + k ( s 

- 1)  > . 

Apply i 

ng 

the  Gauss 

-Green  theorem  and  u 

sing  (4.19), 

(4.16)  i 

we  arrive, 

as 

in 

the 

proof 

of 

Theorem 

4.1,  to 

an 

estimate  of  the  form 

(4.14) 

with 

(4. 

20) 

y(t) 

- 1 

1 (Y-v.F- 

F,  n-n)  ( • 

,t)| 

L2  ( 1 X| <a+  k(s 

-t))’ 

(4. 

21) 

g(t) 

- 1 

1 (b-b, r- 

r) (• ,t)  | 

1 2 

» 

L ( 

| X | <a  + k(s-t) ) 

so  that  Lemma  4.1  yields  (4.12).  The  proof  of  the  theorem  is 
complete . 

From  the  proof  of  Theorems  4.1  and  4.2  it  becomes  clear  that 
in  the  one-dimensional  case,  where  the  internal  energy  is  globally 
convex,  one  may  select  6,  in  (4.8)  or  (4.11"),  arbitrarily  large 
so  that  the  stability  and  uniqueness  results  become  ^obal.  As  a 
matter  of  fact,  in  the  one-dimensional  case  and  um.er  some 
additional  assumptions  (most  notably  genuine  nonlinearity)  , 

DiPerna  [9]  has  established  uniqueness  of  piecewise  smooth  processes 
in  the  class  of  proper  admissible  processes.  DiPerna's  analysis 
relies  on  a number  of  novel  ideas  and  the  problem  of  extending  it 
to  the  multidimensional  case  is  still  open. 
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5.  Stability  of  Smooth  Processes  in  the  Stroni 


Elliptic ity  Region 


The  object  of  the  restrictive  convexity  assumption,  in 
Section  4,  was  to  ensure  that  Il(X,t)  be  pointwise  positive 
definite.  In  reviewing  the  proof  of  Theorem  4.1,  however,  it 

' r 

becomes  clewr  that  definiteness  of  I H(X,t)dX  would  suffice  for 

stability.  We  prove  here  that,  whenever  the  body  is  smooth  and 
the  motion  of  its  boundary  is  prescribed,  stability  is  induced  by 
a mere  strong  ellipticity  condition. 

W'e  will  say  that  a process  (x,n)(X,t)  resides  in  the  strong 
ellipticity  region  if  there  is  a positive  constant  v with  the 
property  that  for  any  vectors  <,  , <;  e Rn,  every  p e R and  all 
(X,t)  in  the  domain  of  the  process 


(5.1) 


aV(F(X,t),n(X,t);X) 


32e*(F(X,t),n(X,t);X) 

+ 3F~5n  5i4<*U 


a2fc*(F(X,t) ,n(X,t) ;X)  - 

♦ ^ = p2  > p(KI2KI2  - p2). 

3n 


Strong  ellipticity  in  nonlinear  elasticity  has  been  studied 
in  connection  to  wave  progagation  (e.g.  [ 1 3 J ) as  well  as  to 
existence  of  equilibrium  configurations  [14].  As  an  assumption 
it  is  weaker  than  convexity  of  the  internal  energy  and,  in  particular, 
it  is  not  incompatible  with  frame  indifference  at  a natural  state. 

For  a discussion  on  the  extent  of  the  strong  ellipticity  region 
in  two-dimensional,  homogeneous,  isotropic,  hyperelastic  bodies, 
see  [15].  See  also  [16]  for  relevant  illuminating  remarks. 


• . 
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I 

i 

f 

We  now  state  the  stability  result: 


Theorem  5.1.  Assume  that  & is  smooth  and  bounded.  Let 
(x,n)(X,t)  be  a smooth  process  defined  on  .0  x [ 0 , t ^ ] , residing 
in  the  strong  ellipticity  region,  with  supply  terms 
(b,r)(X,t)  e L“(.«  x [ 0 , t n ] ) . Then  there  are  positive  constants 
<5,a,3,M>N  with  the  following  property:  If  (x,n)(X,t)  is  any 
admissible  process  defined  on  & x [0,t„],  with  supply  terms 
(b,r)(X,t)  e L 1 ( [0  , tp]  ; L2 ( .0) ) and  such  that 

(5.2)  | F(X,t)  - F(X,t) I + |n(x,t)  - n(X,t)|  < 6 (X,t)  e £ x [0,tQ] , 

(5.3)  x(X,t)  = x(X,t),  X £ 3.0,  t e [0, tQ]  , 
then  we  have,  for  any  s e [0,t^], 

(5.4)  I I (V-V.F-F, n-n) (• ,s) I I 2 r 

L (#) 

< M exp(as+6s2) | | (v-v,F-F,n-n) ( • ,0) | | ? 

, LZ(.0) 

+ N expO»s+fls2)  [ | | (b-F, r-r)  ( • , t)  | , ' 7 dt . 

Jo  ~ ~ L (#) 

t 

From  the  above  proposition  one  draws  the  following  uniqueness 

result : 


► 
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Corollary  5.1.  Let  <#,  (x,n)(X,t)  and  (x,n)(X,t)  he  as  in 
Theorem  5.1.  Assume  that  the  correspond  ing  supply  terms  (b  r)(X,t) 
and  (b,r)(X,t)  coincide  on  x lO.t^l  and  that  both  processes 
originate  from  the  same  state,  i.c., 


(5.5) 


x ( X , 0 ) = x (X , 0) , v(X,0)  = V (X , 0) , X e Q . 


Then  (x,n)(X,t)  and  (x,n)(X,t)  coincide  on  & x [0,t^]. 

For  the  proof  of  Theorem  5.1  we  will  employ  the  following 
Gronwall  type  inequality. 


Lemma  5.1.  Assume  that  the  nonnegative  functions  y(t)  £ L [0,s] 
and  g(t)  e L^[0,s]  satisfy  the  inequality 


(5.6)  y 2 ( o ) < M2y2(0)  + [° [ ( 2y  + 4 fa) y2 (t)  + 2Ng ( t) y ( t) ] dt , o e [0,s], 

J n 


with  P,y,M,N  nonnegative  constants.  Then 


2 2 

(5.7)  y(s)  < M exp(cxs  + Ps  )y(0)  + N exp(<*s+Ps  ) g(t)dt, 

J n 


where  a = y + p/y. 


Proof . We  define  a nonnegative  function  z(o)  by 


(5.8)  z 2 ( o)  = M2y2(0)  + (° ( (2y  + 4po)y2 (t)  + 2Ng (t) y ( t) ] dt , a e (0,s]t 

J n 


and  we  note  that 


(5.!)) 


2 z ( o ) ?.  ( O ) 


2Nfi(o)y(o)  + 4 P 


n 


= (2Y+43o)y'(o)  + 


< (2u  + 4 go)  z“  (o)  + 2Np(ci)z(o). 


y2(t)dt 

0 


Hence 

(5.10)  z(o)  < (<x 2 P o ) z(o  ) +■  N g ( o ) . 

Integrating  the  differential  inequality  (5.10)  under  the  initial 
condition  z(0)  = My(0)  wc  arrive  at  (5.7).  The  proof  is  complete. 

IVe  note  that,  when  ft  = 0 , Lemma  5.1  reduces  to  Lemma  4.1. 


Lemma^  5.2.  Let  (x,n)(X,t)  and  (x,n)(X,t)  he  as  in  Theorem  5.1. 

Then  there  arc  constants  A > 0 and  k with  the  property  that,  for 
any  o e f 0 , 1 0 ] , 


(5.11) 


P(X) 


a2e*(F(X,o),n(X,o);X) 


lFia(X,o)  - Fia(X,o)](L.e(X,o)  - Fjp(X,o)] 


92P*(F(X,o),n(X,o);X) 


(X,o) 


Fia(x,o)Hn(x,o)  - n(x,o)i 


a2t*(i'(x,o),n(x,o);x) 

““17“  “ 


Fn(X.o)  - n(X 


,o)l2jdX 


> 2A 


{ |F(X,o) 

m 


F(X,o)|2  ♦ |n(X,o) 


n(X,o)|2)dX 


k f | x (X , o)  - x(X,o)|2dX. 

1 
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Proof . We  recall  that  (x,n)(X,t)  resides  in  the  strong  ellipticity 
region  and  that  x - x vanishes  on  3 <#  . In  the  case  where 
n(X, t)  - n(X  , t)  vanishes  identically,  (5.11)  reduces  to  the  classical 
Cording  inequality.  In  general,  (5.11)  is  established  by  imitating 
the  proof  of  Cording's  inequality,  as  described,  for  example,  in  [17]. 


Proof  of  Theorem  5. 1 . On  account  of  (2.9), 


* 


1 * 


(5.12)  E 0:,n;X)  - E (P,n;X)  - T.a(I;,ri;X)  (P.a-r.u)  - 0 (P,n;X)  (n-n) 


32t*(F,n;X)  _ _ 32£*(F,n;X) 

3I:ia3ICp  (I'ia'I'ia^F'jB'T'j^  + 3F.(Jo  ja)  (n-n) 

32e*(F.n;X) 


3n 


(n-n)2  + o(l F-F| 2 + | n- n| 2) . 


Combining  (5.11)  with  (5.12)  we  conclude  that  there  is  a positive 
constant  <5  with  the  property  that,  when  (5.2)  is  satisfied, 


(5.13) 


<2 


(p(X)e*(F(X,o),n(X,o);X)  - P(X)E*(F(X,o),n(X,a);X) 


- Tia(5(?’°)»rlCX,o);X)[F.a(X,o)  - F.a(X,o)] 


- P(X)0  (F(X,a),n(X,a);X)[n(X,a)  - r,(X,o)]}dX 


> * f f I F(X,a)  - F(X,o) | 2 + |n(X,a)  - n(X,o)|: 
- k | |x(X,o)  - x(X,o)|2dX. 


}dX 
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t 

( 

i 

Furthermore,  upon  using  Schwarz's  inequality  and  the  Poincare 
inequal i ty , 

I 

(5.14)  f | x (X, o)  - x(X,o)|2dX 


_ _ 2 

| x (X , 0)  - x ( X , 0 ) + fv(X,t)  - v(X,t))dt|  dx 


< 2c  | I*  ( X , 0 ) - F (X , 0)  | 2dX 


/•°f  7 

+ 2o  I v(X,t)  - v(X , t) | dXdt. 

J0  J.#  ~ 

IVe  now  proceed  as  in  the  proo f of  Theorem  4.1.  We  fix 
s £ [0,tQ],  we  integrate  (3.10)  over  # x [0,o],  o £[0,sl,  and  we 
apply  the  Gauss-Green  theorem.  We  use  the  boundary  conditions  (5.3), 
(5.13),  (5.14)  and  the  observation  that  the  right-hand  side  of  (3.10) 
is  of  quadratic  order  in  (v-v ,F-F , n- n , b-b , r- r)  thus  arriving  at 
an  estimate  of  the  form  (5.6)  with 

(5.15)  y (t)  = | | (v-v,F-F,n-n) (• ,t) | | 2 

~ ~ L (.#) 

(5.16)  g(t)  = | | (b-F,r-r) (• ,t) | | 7 

' ~ LZ(^) 


An  application  of  Lemma  5.1  yields  (5.4)  and  completes  the  proof. 


I 
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